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Abstract. We solve numerically the equations of motion for 
the collapse of a shell of baryonic matter, made of galaxies and 
substructure of 10 6 Mq — 1O 9 M0 , falling into the central regions 
of a cluster of galaxies taking account of dynamical friction. We 
calculate the evolution of the expansion parameter, a(t), of the 
perturbation using a coefficient of dynamical friction, r/o, cal- 
culated for a perturbation in which clustering is absent and a 
coefficient rj c i obtained from a clustered one. The effect of the 
dynamical friction is to slow down the collapse (V. Antonuccio- 
Delogu & S. Colafrancesco 1994, hereafter AC) producing an 
observable variation of the parameter of expansion of the shell. 
The effect increases with increasing n and with the increasing 
of clustering. Finally, we show how the collapse time depends 
on ?7o and n c i. 
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1. Introduction 

In the most promising cosmological scenarios, structure for- 
mation in Universe is generated through the growth and col- 
lapse of primeval density perturbations originated from quan- 
tum fluctuations (Guth & Pi 1982; Hawking 1982; Starobinsky 
1982; Bardeen, Bond, Kaiser and Szalay 1986, hereafter BBKS) 
in an inflationary phase of early Universe. Density fluctuations 
originated in the inflationary era are Gaussian distributed and 
their statistics can be expressed entirely in terms of the power 
spectrum of density fluctuations: 



P(k) =< |«5 k | 2 > 
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and pb is the mean background density. 
In biased structure formation theory it is assumed that cosmic 
structures of linear scale R form around the peak of density 
field, <5(x), smoothed on the same scale. Density perturbations 
evolve towards non-linear regime because of gravitational in- 
stability which breaks away from general expansion at a time, 
t m , given by: 



3tt 
32Gp b 



(1 + 5) 



1/2 



(1- 



\3/2 



(4) 



Send offprint requests to: A. Del Popolo, M. Gambera 



where z is the redshift, S is the overdensity within r. When 
S ~ 1 the density perturbation begins to collapse. The col- 
lapse time depends on the characteristic of initial fluctuation 
field such as the average overdensity and on the environment in 
which the perturbation is embedded. This last feature depends 
on the cosmological scenario. 

In this paper, we consider the CDM model (Liddle & Lyth 
1993) based on a scale invariant spectrum of density fluctua- 
tions growing under gravitational instability. The modality of 
the collapse of a perturbation depends on the chosen model. A 
very simple model for accretion of matter in cluster of galax- 
ies was firstly investigated by Gunn & Gott (1972). Gunn & 
Gott's model is an oversimplification of the perturbation evolu- 
tion. In their model they considered only the collapse of spher- 
ical uniform overdense shells of matter surrounding an isolated 
top-hat perturbation embedded within a homogeneous back- 
ground without including the tidal interactions among the shell 
of matter and the external density perturbations and excluding 
the presence of substructure (collapsed objects whose length is 
shorter than that of the main perturbation). This does not 
take into account the previsions of CDM models according to 
which there is an abundant production of substructure during 
the evolution of the Universe. 

In this scenario structure formation proceeds bottom-up 
through gravitational clustering, merging and violent relax- 
ation of small scale substructure (White & Rees 1978) . So the 
matter inside a given region is clumped in a hierarchy of ob- 
jects of various dimensions. The substructure acts as a source 
of stochastic fluctuations in the gravitational field of the pro- 
tostructure inducing dynamical friction (AC) and eventually 
producing a modification of the motion of shells of matter in 
a density perturbation. The galaxies inside a shell of matter 
are subject to the stochastic gravitational field produced by 
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the substructure and their motion undergoes a preferential de- 
celeration in the direction of motion. In particular, dynamical 
friction produces a delaying effect in the collapse of the regions 
with 8 < 10~ 2 inside the perturbation as shown in (AC). The 
final result is that Gunn & Gott's model is inadequate to de- 
scribe the collapse of a spherical perturbation and it requires 
a revision. 

To this aim we used the modified equation of motion of a shell 
of baryonic matter given by AC, in which a frictional force 
that takes into account dynamical friction effects is introduced 
(we study only the component of the dynamical friction due 
to the galaxies belonging to the shell and to the substructure, 
neglecting gravitational bound to the shell) and we solve it by 
numerical integration. 
In this paper we show: 

a) how the expansion parameter, a(t), of the shell is changed 
by the presence of substructure both in the case of an unclus- 
tered and in a clustered system; b) how clustering increases 
the effects of dynamical friction; c) the changes produced by 
dynamical friction on the collapse time of the perturbation and 
the role of clustering in this process. 



2. Modification of Gunn & Gott's equation. 

In a hierarchical structure formation model, the large scale cos- 
mic environment can be represented as a collisionless medium 
made of a hierarchy of density fluctuations whose mass, M, is 
given by the mass function N(M, z), where z is the redshift. In 
these models matter is concentrated in lumps, and the lumps 
into groups and so on. In what follows, we consider a shell 
of matter outside the main body of a perturbation which col- 
lapsed to form a cluster of galaxies. 

The equation of motion of a shell of matter (as previously said 
composed of galaxies and substructure) around a maximum 
of the density field (neglecting tidal interactions and the sub- 
structure), can be expressed in the form: 
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(Peebles 1980, Eq. 19.9), where M is the mass enclosed 
with the proper radius r(t). Using Gunn & Gott's notation the 
proper radius can be written as: 



r(n,t) = a(r x ,t)n 
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where r% is the initial radius and a(ri,t) is the expansion 
parameter of the shell. At the initial time t, the initial condition 
is given by 



a(n,ti) = 1 



(7) 



In the presence of substructure it is necessary to modify the 
equation of motion, Eq. (g), because the graininess of mass 
distribution in the system induces dynamical friction that at 
last introduces a frictional force term in Eq. (|^). 
In a material system, the gravitational field can be decomposed 
into an average field, Fo(r), generated from the smoothed out 
distribution of mass, and a stochastic component, F stoc h(r), 
generated from the fluctuations in number of the field particles. 
The stochastic component of the gravitational field is specified 
assigning a probability density, W(F), (Chandrasekhar & von 



Neumann 1942). In an infinite homogeneous unclustered sys- 
tem W(F) has been given by Holtsmark distribution (Chan- 
drasekhar & von Neumann 1942) while in inhomogeneous and 
clustered systems W(F) has been given by Kandrup (1980) 
and Antonuccio-Delogu & Barandela (1992) respectively. The 
stochastic force, F stoc h, in a self-gravitating system modifies 
the motion of particles as it is done by a frictional force. In 
fact a particle moving faster than its neighbours produces a 
deflection of their orbits in such a way that the average den- 
sity is greater in the direction opposite that of traveling causing 
a slowing down in its motion. 

Following Kandrup (1980) in the hypothesis that there are no 
correlations among random force and its derivatives the fric- 
tional force is given by: 



F = -77V = - 
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where r\ is the coefficient of dynamical friction, T(F) the 
average duration of a random force impulse, < v 2 > the charac- 
teristic speed of a field particle having a distance r ~ {^jr-) 1 ^ 2 
from a test particle (galaxy) . This formula permits to calculate 
the frictional force for inhomogeneous systems when W(F) is 
given. If the field particles are distributed homogeneously the 
dynamical friction force is given by: 
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(Kandrup 1980), where m a and n a are respectively the av- 
erage mass and density of the field particles. Using the virial 
theorem we can written the dynamical friction force as follows: 
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where N is the total number of field particles. In this last 
equation it is assumed that the field particles generating the 
stochastic field are virialized. This is justified by the previrial- 
ization hypothesis (Davis & Peebles 1977). 
To calculate the dynamical evolution of the galactic component 
of the cluster it is necessary to calculate the number and aver- 
age mass of the field particles generating the stochastic field. 
The protocluster, before the ultimate collapse at z ~ 0.02, is 
made of substructure having masses ranging from 10 6 — 10 9 A/q 
and of galaxies. We suppose that the stochastic gravitational 
field is generated from that portion of substructure having a 
central height v larger than a critical threshold u c . This latter 
quantity can be calculated (following AC) using the condition 
that the peak radius, r v k(v > u c ), is much less than the aver- 
age peak separation l av = n a (y > v^)~ x (see AC and BBKS). 
The condition r v u{v > v c ) < 0.1n a (^ > fc)^ 1 ^ 3 ensures that 
the peaks of substructure are point like. Using the radius for a 
peak (AC Eq. 13) that is: 
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where 7, R* are parameters related to moments of the 
power spectrum (BBKS Eq. 4.6A), we obtain a value of 
v c = 1.3 and then we have n a (y > v c ) = 50.7Mpc -3 (7 = 0.4, 
R* = 50Kpc) and m a is given by: 
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(in accordance with the result of AC), where m p k is given 
in Peacock & Heavens (1990) and N p k is the average number 
density of peak (BBKS Eq. 4.4). Now we study the evolution 
of a shell of matter, made of substructure and galaxies, as mod- 
ified by dynamical friction. 

As previously told, we suppose that the substructure has a 
Maxwellian distribution of velocity and that the galaxy moves 
in the substructure background. The modified equation of mo- 
tion for each galaxy of a shell of matter can be written in the 
form: 



dt 2 



GM 
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(Langevin 1908, Kandrup 1980, Saslaw 1985, Kashlinsky 
1986, AC) In the spherical collapse model, it is supposed that 
the infall is radial for every galaxy in the protocluster. Eq. ( ^) 
describes the motion of each galaxies, the only parameter that 
describes the collapse is the radial distance, r, or the parameter 
a(t) previously defined. Remembering that the average density 
is given by: 



P(ri,t) 
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and that mass conservation requires: 
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equation ( |l^ ) in terms of a(r%,t) can be written as: 
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where p C i is the background density at a time tj and Si 
is the overdensity within n. Using the parameter r = 
where T c o/2 is the collapse time in absence of dynamical fric- 
tion (Gunn & Gott Eq. 16), Eq. ( [t6| ) may be written in the 
form: 



d 2 a 
dr 2 
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Equation ( |17[ ) is obtained remembering that the probabil- 
ity density W(F) depends linearly on the correlation function 
(AC) and so it is possible to decompose the coefficient of dy- 
namical friction, 77, as follows: 



V = rj + T]d 



(18) 



where 770 is the coefficient of dynamical friction of an un- 
clustered distribution of field particles, while r\ c i takes into ac- 
count clustering. 

3. Results and discussions. 

The time evolution of the expansion parameter, a(r), can be 
obtained solving Eq. (p^). 

There are two ways to solve equation (|l7|). The first one is to 



look for an asymptotic expansion, the second one is numeri- 
cal. In AC the authors gave only an expression for the collapse 
time, T c , taking into account dynamical friction. We solved Eq. 
( |l7j ) numerically using a Runge-Kutta integrator of 4 th order 
and we studied the motion of a shell of matter of low density, 
S — 0.01, typical of a perturbation present in the outskirts of 
a cluster of galaxies. We chose the initial conditions remem- 
bering that at the maximum of expansion the initial velocity 
is zero. We show in fig. 1 the results of our calculation when 
the dynamical friction is present but considering only the case 
of an unclustered distribution ( 7/0 7^ and rj c i — 0). 
In this case the expansion parameter of a shell, a(r), is plotted 
versus r, for different values of 770, the coefficient of dynamical 
friction for an unclustered system (stochastic force generators 
randomly distributed). 



Fig. 1. Temporal evolution of the expansion parameter of a shell 
of matter made of galaxies and unclustered substructure. The solid 
line is a(r) when dynamical friction is absent while the dotted line 
is a(r) when it is taken into account ( 1)0 / and tj c ; = 0). The 
dashed line shows the effect of dynamical friction when 77^ > 770. We 
assume a cluster radius of R c i = 5h~ Mpc, a central overdensity 
<5 = 0.01 and a total number of peaks of substructure Ntot = 10 3 . 

The effect of dynamical friction is to make the decrease of 
a(r) less steep and consequently to slow down the collapse of 
the shell in agreement with the analytic calculation performed 
by AC. The evolution of the parameter of expansion, o(t), 
is more and more modified as 770 increases, in comparison to 
Gunn & Gott's model, as one can see from figure 1. 
In Fig. 2 we show the same parameter, o(t), for different values 
of ?7 c z, the coefficient of dynamical friction of a clustered sys- 
tem (clustering of substructure is required, for example, by the 
minihalo model (Mo et al. 1993)). Clustering produces a sen- 



Fig. 2. Temporal evolution of the expansion parameter of a shell of 
matter made of galaxies and clustered substructure. The solid line 
is o,(t) when dynamical friction is absent while the dotted line is 
the same when it is taken into account (r; c ; ^ 0). The dashed line 
shows the effect of dynamical friction when 7/', > 77 c j. We assume 
a cluster radius of R c i = 5?t — 1 Mpc, a central overdensity 8 = 0.01, 
and a total number of peaks of substructure Ntot = 10 3 



sible slowing down of the shell collapse with respect to unclus- 
tered systems (?7 C ; = 0). This result is perfectly understandable 
since the introduction of positive clustering produces a greater 
probability for a test particle to be scattered during its motion, 
increasing the role of dynamical friction. 

The caracteristic values for 77 are between < 77 < 3, they 
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depend strongly on the value ol many parameters like v c and 
others (Gambera et al. 1996, in preparation). For a typical 
protocluster configuration like that ol Coma, that can be con- 
sidered almost as a lower limit having an initial overdensity, 
at 2 ~ 10 3 , given by 6 ~ 8 • 10 -3 and a total mass estimated 
in M ~ 1O 15 M0 (Hughes 1989), supposing that small scale 
substructure is made of objects of M ~ 1O 9 M0, we obtain 
r\ ~ 0.2. 

Finally, we computed the time of collapse of a shell, T c , versus 
both r/o and of r] c i ■ To this aim we solved Eq. ( |l7| ) varying r/o 
and for each value of this parameter we obtained the zero of 
a(r), that gives us T c . We repeated the same procedure for r\ c i. 
In Fig. 3 we show the results of the computation. 
It is interesting to observe that T c is very sensitive to the 
changes of the coefficent of dynamical friction for a clustered 
system n c i and increases with increasing rj c i while it is less 
sensitive to the changes of the coefficent of dynamical friction 
for an unclustered system. Dynamical friction changes the dy- 
namical evolution of low density perturbations (5 ~ 0.01) and 
clustering contributes to increase this effect. 



Fig. 3. The variation of time of collapse r 



with r\. The 



solid line is the time of collapse of an unclustered system in which 
dynamical friction is present while the dotted line is the same for 
a clustered system and the dashed line is the time of collapse for a 
system in which dynamical friction is absent 
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4. Conclusions 

In this paper we have discussed the dynamics of the infall of a 
shell of baryonic matter of a spherical perturbation, when the 
effects of dynamical friction are taken into account, using the 
time evolution of the parameter of expansion, o(r). In particu- 
lar, we showed that the parameter a(j) of a shell of matter of a 
protocluster decreases less steeply with increasing r/o and that 
clustering produces a further increase in this effect. This effect, 
then, could be the cause because many clusters of galaxies are 
not yet relaxed. Finally, we showed how the collapse time, T c , 
varies in presence or absence of dynamical friction (both in the 
case in which r\ c i = and when n c i 7^ 0). 

The collapse time, T c , of an infalling shell increases with in- 
creasing values of 770 (see Fig. 3) . Clustering enlarges the change 
in the dynamical evolution of the perturbation produced by 
dynamical friction. An interesting point on which we are in 
progress (Gambera et al. 1996, in preparation) is the determi- 
nation of how the growth of the collapse time depends on v c 
and on others parameters. Besides, we want to find an analytic 
relation that links r with the parameters on which it dependes. 
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